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^ 1 Abstract 

H . ^ 
Qh' The gluon density in nuclei, G (x), is poorly constrained at all x from the 

current DIS and Drell-Yan data. In this paper, we point out that J/ip production 
measured in proton-nucleus collisions at y / s^" = 38.8 GeV by the E866 collab- 
oration puts stringent constraints on the x-dependence of the G w /G Be ratio in 
the x = 2 10- 2 -10 _1 range. The E866 data suggest a rather mild x-dependence 
of G w /G Bc , and consequently tend to favour DS and HKM sets, rather than 
the DSg, EKS, and EPS parametrizations which exhibit a large shadowing and 
anti-shadowing. 
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A lot of effort has been devoted over the past ten years in order to constrain the par- 
ton distributions in nuclei. The first set of nuclear parton densities (nPDFs) at leading 
order (LO) was given by Eskola, Kolhinen, Ruuskanen and Salgado (EKS) [1]. This 
parametrization used deep inelastic scattering (DIS) data to determine the valence- 
quark densities, and Drell-Yan production in proton-nucleus collisions to constrain the 
sea. As first done in [2], the gluon sector in EKS is probed indirectly through the 
Q 2 -dependence of the structure function ratios measured by the NMC collaboration. 
However, the rather large error bars of those data could not allow for a precise deter- 
mination of the gluon nuclear distributions. The EKS set has been widely used ever 
since to predict the nuclear dependence of a large variety of hard observables: Drell- 
Yan, electroweak bosons, large-p ± hadrons, heavy quark and heavy-quarkonium, jets, 
and prompt photons. A few years later, Hirai, Kumano, Miyama and Nagai (HKM) 
proposed an alternative distribution, based on an explicit x 2 analysis of the data [3] 
The kinematic dependence of the EKS and HKM sets is shown to be rather different, 
as discussed e.g. in [4], reflecting the relative lack of constraints given by the available 
data. The first global fit analysis of nPDFs at next-to-leading order (NLO) accuracy 
has then been performed by De Florian and Sassot (DS) [5]. Finally, an analysis taking 
into account non-linear corrections in the QCD evolution at small x is also proposed 
in [6]. 
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In order to further constrain the gluon nuclear density - whose knowledge is essential 
to perform reliable predictions in proton-nucleus and nucleus-nucleus collisions at the 
LHC -, Eskola, Paukkanen, and Salgado (EPS) recently included in their analysis [7] 
the forward hadron production data measured in d-Au collisions by the BRAHMS 
collaboration at RHIC [8]. The rather large suppression in d-Au with respect to 
p-p collisions seen experimentally leads to a dramatic effect in the ratio R G (x) = 
G A (x)/G p (x) <C 1 at small x = O (10 -3 ) and, following momentum sum rules, to a 
significant anti-shadowing, R (x) ~ 1.4 at x — (9(10 _1 ) and at rather low scales. As 
a consequence, the x-dependence of G A /G p (x) in EPS proves much steeper than the 
one predicted in any other nPDF sets. 

In this paper, we investigate J / ip production in proton- nucleus collisions as an inter- 
esting channel to constrain the gluon nuclear density. Indeed, assuming that it follows 
that of open charm, J /if; production proceeds via gluon fusion and quark- ant iquark 
annihilation. However, as long as J /if; is produced at small longitudinal momentum 
fraction, x-p <C 1, the gluon fusion channel dominates over the qq annihilation process. 
At leading order, the differential cross section is therefore simply proportional to the 
product of gluon densities: 

' 1 " (p A - J /if; X) oc G p (x 1 , Q 2 ) G A (x 2 , Q 2 ) 5{ Xl x 2 s - m 2 ), (1) 



where 



*v = \ (^+4m^/ a ±x P ) (2) 

are the projectile and target-parton momentum-fractions (yfs being the centre-of-mass 
energy of the hadronic collision), and the typical scale Q at which the gluon densities 

should be evaluated is given by the J/ if; mass scale, Q = O (rrij^. In this kinematics, 

the nuclear production ratio, 

R A/P (x 2 ) = j ^(p + A^J/iP + X)/^-(p + p^J/il; + X), (3) 

therefore reduces to R A/p (x 2 ) ~ R G (x 2 ), which we would like precisely to constrain. 
This relationship was first used in [9] to extract the gluon nuclear density from J /if) 
data in 7r-A and p-A collisions measured by the NA3 and E772 collaborations. 

However, J /if) is a pretty weakly bound state which may be sensitive to inelastic 
rescattering processes in large nuclei, which will spoil the above simple relationship 

A 

between R and R Q . Assuming the factorization between the cc production and the 
rescattering process, the nuclear production ratio can be written as: 

R A/P M~ R g( x *) x ^abs(A,a WN ), (4) 

where 5 , a b s (A, o" 7/l ^ N ) denotes the probability for no interaction (or "survival probabil- 
ity") of the J /if) meson with the target nucleus. It depends of course on both the 
atomic mass number A of the nucleus and the J/ifi-N inelastic cross section, cr N , 
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and is given in a Glauber model by [10] 



&»(A, gjwl ) = (A _ i — /db (l - «-<-""> W . w ) , ( 5 ) 

where T A (b) is the so-called nuclear thickness function. Because of the present lack of 
constraints on cr N (and therefore on S^), the hope to determine R Q (x) from J/ip 
suppression data in p-A collisions is rather limited. It is the reason why this channel 
has never been considered in the global fit analyses of nPDFs. As we shall see later, 
significant constraints can nevertheless be achieved by fitting the J / ip suppression data, 
letting the J / ip inelastic cross section as a free parameter. 

The present analysis relies in particular on the measurements performed by the 
E866 fixed-target experiment in p-Be and p-W collisions at ^/s^ = 38.8 GeV. The 
E866 collaboration reported on the nuclear production ratio [11]: 

R w/Bc (x F ) = ^ ^(p + W -> J/tf + X) /^-(p + Be - + X) (6) 

where A Bc = 9 and A w = 183 denote the atomic mass numbers. The xp range covered 
by these data is fairly large, —0.1 < xp < 0.9. At large xp, R w/Bc decreases very 
rapidly down to 0.3 in the largest x-p bin. The reason for such a large suppression is not 
settled yet. In any case, the J/ip inelastic interaction together with the modifications 
of parton densities in nuclei - the only effects assumed in the present work - could 
explain only partly this large-^F suppression [12]. Furthermore, the lack of x 2 -scaling 
when comparing large-xp E866 data points with NA3 [13] and PHENIX [14] data 
indicates that other effects, beyond inelastic rescattering and nPDFs corrections, are 
responsible for the suppression reported at large xp. Therefore, the present analysis will 
be restricted to J/ijj production around mid-rapidity, \xp\ < 0.25, which corresponds 
to the range x 2 = 2 10 _2 -10 _1 in the target momentum-fraction. 

Even though Eq. (jlj) should be a good approximation for the J ftp suppression in 
proton-nucleus collisions, the calculation is performed in the colour evaporation modefj], 
including both the gg and qq channels for the production of cc pairs with invariant 
masses ttLj/^ < m < 2mo (see Ref. [16] for more detail on the present calculation). 
The W/Be nuclear production ratio measured by E866 has been fitted using the various 
nPDFs currently available, DS [5], DSg [5], EKS [1], EPS [7], HKM [3], keeping a J/iin 
as a free parameteJl. It is not the goal to discuss the actual value of <J J/iiN obtained, 
we refer the reader to [16] for a comprehensive discussion on this point. Rather, the 
respective agreement between the fits and the E866 measurements using each nPDF set 
is investigated. Moreover, because of the above mentioned nuclear absorption which 

lr Tliis model describes heavy-quarkonium production near threshold and reproduces fairly well 
rapidity and low-p ± spectra in hadronic collisions, although it is known to face some difficulties (see 
e.g. [15] for a critical discussion). In the present context, this approach is only used to estimate the 
small corrections to Eq. ([4]) because of the quark annihilation process; hence our results do not depend 
much on which specific model is assumed. 

2 For the time being Cj,^ n is taken to be independent of x 2 , although we shall discuss in the 
following the implications of a possible cc 2 -variation of the cross section on our results. 
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strength is not precisely determined, these data will rather probe the ^-dependence 
of R Q (x 2 ) rather than its absolute value. 
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Figure 1: The W-over-Be J/ip production ratio as a function of x 2 measured by the 
E866 collaboration, in comparison with the best fits using the proton, DS, HKM (left) 
and DSg, EKS, EPS (right) parton densities. The overall 3% systematic errors (i.e. 
independent of not shown. 

As can be seen in Fig. [TJ the E866 data points show a remarkably flat behaviour as 
a function of x 2 . Therefore, following Eq. (@J, the gluon nuclear density ratio R G (x) at 
the J/ip mass scale seems to have a pretty limited variation, say < 5% from x = 2 10~ 2 
to x = 10 _1 . Let us compare this observation with the predictions given by the 
different nPDF sets. The trend reported by the E866 data is well reproduced by the 
DS and HKM sets on the entire x 2 -domain. These two parametrizations indeed do not 
predict a strong ant i- shadowing at x 2 = O(10 _1 ) nor a large shadowing at small x 2 . 
Moreover, note that the E866 measurements are also consistent with no modifications 
of the gluon distributions in nuclei, R G (x) = 1 (labelled proton), or at least with an 
^-independent ratio, R G (x) = R G , because of the uncertainty of the normalization due 
to the nuclear absorption. The good agreement between the data and these theoretical 
calculations is also indicated in Table [JJ where the xV n df for each calculation is given 
and turns out to be O (1). On the contrary, the nPDFs sets which predict a strong 
dependence^] of R G on x, like DSg, EKS, and EPS, seem to be somehow disfavoured by 
these measurements. In the fits shown in Fig. [1] (right), the large (EKS, DSg) and very 
large (EPS) variation assumed in those sets turns out to be in striking contradiction 
with the data, as also indicated by the large values x 2 /ndf (respectively 10.1, 12.0, and 
35.9) given in Table [TJ It is remarkable that the strong shadowing present in EPS in 
order to reproduce the large-rapidity hadron production data measured by BRAHMS 
at RHIC appears to be too large to reproduce the J/ip data in the x 2 = 10~ 2 — 10 _1 
range. It may indicate that it is not possible to fit simultaneously the forward RHIC 

3 More precisely a much stronger dependence of with respect to that of R® c - 
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results and lower energy (larger x) data in terms of a universal modification of parton 
densities in nuclei. 

We investigated as well the constraints given by the NA3 measurements in p-A 
collisions at lower beam energy, A /s^" = 19.4 GeV [13]. However, the limited x 2 ~ 
7 10 _2 -10 _1 coverage (restricting again to not too large xp) as well as the larger error 
bars do not allow for differentiating the various nPDFs, all predictions leading to 
X 2 /ndf=0(l). 



nPDF set 


proton 


DS 


DSg 


EKS 


EPS 


HKM 


xViidf 


0.5 


0.6 


12.0 


10.1 


35.9 


1.2 



Table 1: xV n df values between the best fits and the E866 data using the proton, DS, 
DSg, EKS, EPS, and HKM nuclear parton densities. 

As already discussed, one should however keep in mind that the nuclear absorption 
is an important process at work in the J/ip suppression. In particular, one could 
argue that a fast variation of the nuclear absorption as a function of x 2 may actually 
compensate the apparently too strong x-dependence of R G (x) assumed in DSg, EKS, 
and EPS. In perturbative QCD, cr //</)N depends slightly on the incident J/^-nucleon 

center of mass energy, cr J/iin oc \ JsJ^\ , where the exponent A ~ 0.25 is related 

to the behaviour of the small- x gluon distribution in the nucleon, xG(x) ~ x _A [17]. 
Since y/s J/iiN ~ m JW /^, the cross section is expected to scale like x~ A . Expanding 
the exponential in (j^J) leads to S a b B (x 2 ) ~ 1 — #/x A which increases with x 2 . The 
energy dependence of o" J/v , N would actually worsen the agreement between the fits and 

the data. The only possibility to reconcile the rapid x-dependence of R with the flat 
behaviour of the E866 data would be to appeal to dramatic formation time effects: when 
x 2 is large, the J/ip state is formed rapidly (recall that the Lorentz boost 7 oc l/x 2 ) 
and could experience a stronger interaction with the nucleus than a more compact cc, 
at small x 2 , propagating through the nuclear medium. The competition between such 
formation time effects together with the strong x-dependence of R (x) may then result 
into an accidentally flat ratio R w/Bc - This explanation would be really fortuitous and 
appears therefore rather unlikely. We do not really see either how NLO corrections 
could alter the present conclusions. Because of the extra gluon radiation, the gluon 
densities could be probed at slightly larger values of x 2 than what is estimated in 
Eq. ([2]) at LO. This would nevertheless not improve the agreement between the data 
and the DSg, EKS, and EPS fits. 

To summarize, we discuss in this paper how J/ip suppression measurements per- 
formed in proton-nucleus collisions {y/s~^ = 38.8 GeV) by the E866 experiment [11] 
gives interesting constraints on the x-dependence of the gluon distribution in nuclei. 
The flat behaviour observed experimentally suggests a small variation of the G A /G p (x) 
ratio in the x = 10~ 2 -10 _1 range at the J/ip mass scale, which is consistent with DS 
and HKM parametrizations. On the contrary, the large shadowing and antishadowing 
effects (and from this the fast variation of G A /G p (x)) assumed e.g. in DSg, EKS, and 
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EPS seem to be disfavoured, unless fortuitous effects due to the interaction of dynam- 
ical cc pair in the nuclear medium are present. In that sense, comparing the present 
E866 J / ip data with the isolated prompt photon measurements in d-Au collisions soon 
to be taken at RHIC in the p ± = 2-5 GeV range (probing the gluon nPDF at roughly 
to x 2 = 4 lO^-lCT 1 [18]) will be particularly interesting. 
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